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Abstract 

We consider a phase field model for the flow of two partly miscible 
incompressible, viscous fluids of Non-Newtonian (power law) type. In 
the model it is assumed that the densities of the fluids are equal. We 
prove existence of weak solutions for general initial data and arbitrarily 
large times with the aid of a parabolic Lipschitz truncation method, 
which preserves solenoidal velocity fields and was recently developed 
by Breit, Diening, and Schwarzacher. 
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1 Introduction 

We consider the flow of two macroscopically immiscible, incompressible Non- 
Newtonian fluids. In contrast to classical sharp interface models, a partial 
mixing of the fluids is taken into account, which leads to a so-called diffuse 
interface model. This has the advantage that flows beyond the occurrence 
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of topological singularities e.g. due to droplet collision or pinch-off can be 
described. More precisely we consider 



p^tv + pv ■ Vv - div S(c, Dv) + Vp = div(Vc (g) Vc), (1.1) 

divv = 0, (1.2) 

(9fC + V ■ Vc = mAp, (1.3) 

^ = £-V(c) - eAc (1.4) 

in Qt = Q X (0,T), where Q C R"', n > 2, is a bounded domain and 
T G (0, cxd). Here v is the mean velocity, Dv = |(Vv + Vv"^), p is the 
pressure, c is an order parameter related to the concentration of the fluids 
e.g. the concentration difference or the concentration of one component, and 
p is the density of the fluids, which is assumed to be constant. Moreover, 
S(c, Dv) is the viscous part of the stress tensor of the mixture to be specified 
below, e > is a (small) parameter, which is related to the "thickness" 
of the interfacial region, $ : M — )• M is a homogeneous free energy density 
and = $' and p is the chemical potential. Capillary forces due to surface 
tension are modeled by an extra contribution eVc (g) Vc := £:Vc(Vc)^ in the 
stress tensor leading to the term on the right-hand side of (11.11) . Moreover, 
we note that in the modeling diffusion of the fluid components is taken into 
account. Therefore mA/i is appearing in (11. 3p . where m > is a constant 
mobility coefficient. 

We close the system by adding the boundary and initial conditions 

v|9n = ondQx{0,T), (1.5) 

n ■ Vclan = n ■ Vp|af7 = on dil x (0, oo), (1.6) 

(v,c)|i=o = (vo,co) in 1^. (1.7) 

Here n denotes the exterior normal at dQ. We note that (11.11) can be replaced 
by 

pdtv + pv ■ Vv - div S(c, Dv) + = pVc (1.8) 
with g = p + f I Vcp + £"^$(0) since 

;xVc = V (|l Vc|2 + £'^$(c)) - 6 div(Vc ® Vc). (1.9) 

In the case of Newtonian fluids, i.e., S(c, Dv) = z/(c)Dv for some pos- 
itive viscosity coefficient z^(c), the model was first discussed by Hohenberg 
and Halperin [13]. Later it was derived in the frame work of rational con- 
tinuum mechanics by Gurtin, Polignone, Vinals [12]. The latter derivation 
can be easily modified to include a suitable non-Newtonian behavior of the 
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fluids. If e.g. S(c, Dv) is cliosen sucli tliat S(c, Dv) : Dv > 0, tlie local dissi- 
pation inequality, which yields thermodynamical consistency, remains valid. 
For results on existence of weak and strong solutions in the case of Newto- 
nian fluids we refer to Starovoitov [19], Boyer [6], and A. |3]. First analytic 
results for the system fll.ll) - fll.4p for Non-Newtonian fluids of power-law type 
were obtained by Kim, Consiglieri, and Rodrigues [2]. The authors proved 
existence of weak solutions if g > d = 2,3, where q is the power de- 

scribing the growth of the stress tensor with respect to Dv. For this range of 
q monotone operator techniques can be applied. Moreover, in the case c? = 3 
and 2 < g < y the authors prove existence of weak solutions. Grasselli and 
Prazak [10] discussed the longtime behavior of solutions of fll.ip - fll.4p in the 
case q > = 2,3. 

The goal of this article is to extend the existence result to lower values of 
q in order to include the physically important case of shear thinning flows. 
In the case of a single fluid existence of weak solutions for power-law type 
fluids was proved for the case q > d > 2, by D., Ruzicka, and Wolf [7]. 
The proof is based on a parabolic Lipschitz truncation method and a careful 
decomposition of the pressure, which is needed since the Lipschitz truncation 
used does not preserve the divergence freeness of a velocity field. Recently a 
parabolic Lipschitz truncation method, which keeps divergence free velocity 
fields divergence free, was developed by Breit, D. and Schwarzacher [9]. In the 
present article we will use this method in order to prove existence of weak 
solutions to fll.ip - fll.7p if S(c, Dv) is of power law type with an exponent 
q > -j^. Precise assumptions are made in the following. 

For simplicity we assume that e = p = 1. But all results are true for 
general (fixed) e, p > 0. Moreover, we assume: 

Assumption 1.1 Let C M'^, d = 2,3, be a bounded domain with C^- 
boundary and let $ G C([a, 6]) fl C^((a, h)) be such that (p = ^' satisfies 

lim0(s) = — oo, lim0(s) = oo, (p'{s) > —a 

for some a eR. Let m > and let S: [a, b] x R'^^'^ — )■ R'^^'^ be such that 

|S(c,M)| < C(|sym(M)|''-i + l) (1.10) 
|S(ci,M)-S(c2,M)| < C|ci-C2|(|sym(M)|''-i + l) (1.11) 
S(c,M):M > /t|sym(M)|''-Ci (1.12) 

for all M G M"'^'^, c, Ci,C2 G [a,b], and some C,Ci,K > 0, q E ij^^oo). 
Moreover, we assume that S(c, ■): Rgy^ "^tym ^■^ strictly monotone for 
every c G [a, b] . 
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For the following we denote 



Emix{c) = I ——dx+ / Lp{c)dx. 
Jn ^ Ju 



Let V G L<i{{),T-WlQ{SlY) H L~(0, T; c G L°°{0,T; H\Q)) n 

L\0,T;H\Q)) with (/?(c) G ^^(fi x (0,r)), and /x G L\0,T; H^Q)), where 
< T < oo. Then (v, c, /i) is a weak solution of the system fll.ip -( ITT|) if for 
any ip G C°°{QtY with div^s = and supp((/)) CC Vt x [0,T) the following 
holds: 

-/ Y -dtcfi d{x,t) - V® V : Dc^ + / S(c, Dv) : Dc^ t) 

^Qt JQt JQt 

= e Vc^Vc-.Bcp d{x,t) + / vq ■ (^(0) rfx (1.13) 
and for every 4j e C^{Q x [0, T]) with supp CC H x [0, T) 



n 



cdtip d{x,t) — I 001^(0) dx+ I {v ■ Vc)^jj d{x,t) 

Qt 

-m I Vfi-Vip d{x,t), (1.14) 



/i = 0(c) — Ac, 
n ■ Vc\en = 



holds. 



THEOREM 1.2 Let Assumption fO feo/d trae and /et < T < cx). T/ien 
/or any vq G i^^(^^) and cq G if"'^(il) wt/i Co(a;) G [a,b] almost every- 
where there exists a weak solution v G L''{0,T;W^Q{n)) n L°°(0, T; (fi)), 
c G L°°(0,T; //!((})) n L2(0,r;i/2(^)) y;^^/^ ^(c)' e /.^(fi x (0,r)), and 
fi G L^(0, T; iJ^(fi)) m i/ie sense above. 

The structure of the article is as follows. In Section [2] we summarize 
some preliminary results needed in the following analysis. In Section [3] prove 
existence of solutions of an approximate system to fll.ip - fll.7l) . where the 
convective terms v ■ Vv, v ■ Vc, and the capillary term div(Vc ® Vc) are 
smoothed in a suitable way. Finally, in Section H] the existence of weak 
solutions is proved by passing to the limit in the approximate system with 
the aid of a solenoidal parabolic Lipschitz truncation method. 
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2 Preliminaries 



Throughout the paper the usual Lebesgue spaces with respect to the Lebesgue 
measure are denoted by L^(M), 1 < p < oo, for some measurable M C M^. 
Moreover, U'i^M] X) denotes its Banachspace- valued variant and L^(0, T; X) = 
L^((0,T);X). Furthermore, /: [0,oo) ^ X is in Lf^^([0, oo); X) if / is 
strongly measurable and /|[o,t] £ Lp{0,T] X) for every < T < oo and 

Ll,([0,oo);X) := {/ G ^^([0, oo); X) : ||/|U.^^^([o,oo);X) < oo} 



L5;i^^([0,oo);X) := sup ||/||lp(t,T+1;X)- 

r>o 

The standard L^-Sobolev space is denoted by W^{Q). W^q^Q) is the closure 
oiC^{n) in W^{n) and i7™(n) = W^{n),H^{n) = W^oi^). Furthermore 
we use the notation -^^^o)(^) = {/ ^ -L^(fi) : J^f{x)dx = 0}, iJ(Q)(n) = 
H\n) n Lf^^in), and //(o}(^^) := ^(o)(fi)'- Finally L2(^]) is the closure of 
divergence free C^(fi)-vector fields in L^(f2)'^. 

We recall some results on the Cahn-Hilliard equation with convection 
term: 

dtC + Y ■Vc = mA/i inf^x(0,oo), (2.1) 

fi = (f){c)-Ac infix(0, oo), (2.2) 
n- Vc|an = n- VyuloQ = on (9fi x (0, oo), (2.3) 

c\t=o = Co inQ (2.4) 

for given cq with E^i^{co) < oo and v G L°°(0, oo; Ll{n))nL^{0, oo; ii'o(^)'')- 
Here = $' and $ is as in Assumption 11.11 In the following (12.1 p together 
with (12. 3p will be understood in the following weak form 

{dtc{t) , (f) fj-i ^1 + / v{x,t) ■ Vc{x,t)(p{x) dx = —m / V fi{x,t) ■ \/(p{x) dx 

(0)' (0) 

for all ip G H^Q-j{Q) and almost every t G (0,T), where (■, ■)x',x denotes the 
duality product. Moreover, Qt := 0, x (0,t), Q = Q x (0, oo). 

THEOREM 2.1 Lei v G L\0, oo; i/o^(^^)'^) n L°°(0, oo; Ll{n)). Then for 
every Cq G H^qs^{Q) with Emix{co) < i/iere a unique solution c G 
5(:7([0,oo);f^J))(^])) of [2l\)-[K4\) with dtc G ^^(0, oo; ff(oJ(fi)) an(i /i G 
-^uioclp! oo); iL^(fi)). This solution satisfies 

Emix{c{t))+ I m\Vfi\'^ d{x,T) = E^i^{co) - v ■ fiVcd{x,T) (2.5) 

JQt JQt 
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for all t E [0, oo) and 

|2 , ||;5^||2 , ||V7,,I|2 

(0) 



< C (^ErmAco) + Mh^Q)) (2.6) 

|2 



(2.7) 

where r = 6 if d = 3 and 1 < r < oo is arbitrary if d = 2. Here C, Cr are 
independent o/v,co. Moreover, for every R> the solution 

ceY:= LlM oo); W^m n HlM oo); H^^^m 

depends continuously on 

(Co,v) e X := H\n)xLl^{[{),(x)-Ll{D)) with Emix{Co) + \ML^(^0,oo;m) < R 

with respect to the weak topology on Y and the strong topology on X . 

The theorem is proved in p], Theorem 6] in the case m = 1. The case m > 
can be reduced to the case m = 1 by a simple scahng in time and a scahng 
of the homogeneous free energy density /. 

We need the following theorem, which is a summary of Theorem 2.14 and 
Corollary 2.15 of [9]. 



Theorem 2.2 Let Iq be an open time interval, let Bq be a ball in M.'^, and let 
Qo := Jo X Bq. Let g, a G (1, oo) with q,q' > a > 1, where q' = Let ( G 

C^i^Qo) with xiQo < C < XiQo- "™ "'^^ safe/y dtWm = - div 
in the sense of distributions V'^^^{Qq), where V^n = G C^^QqY : div ip = 
0}. Assume that is a weak null sequence in L'^{Io]W^''^{Bo)) , a strong 
null sequence in L"{Qq) and bounded in L°°{Iq,L"{Bq)). Further assume 
that Gm = Gi^m + G2,m such that Gi^m is a weak null sequence in (Qo) 
and G2,m converges strongly to zero in L'^{Qq). Then there exist a double 
sequence of open sets Om,k, k,m E N, with limsup^_j.oo |Om,fc| < c2~^2~^^'' 
for all k eN such that for every K G (|Qo) 



lim sup 



/ ((Gi,„ + K) : Vu„)CXoC d{x,t) 



< c2-'=/^. 



The proof of Theorem 12.21 is based on a solenoidal Lipschitz truncation. In 
particular, the set Om,k are level sets of suitable maximal operators defined 
by and Gm- It is a refinement of the parabolic Lipschitz truncation of [TT] 
and [7]. The advantage of Theorem 12. 2l is that the pressure can be completely 
avoided by using a solenoidal Lipschitz truncation. See |8j for a solenoidal 
Lipschitz truncation in the stationary case. 
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3 Approximate System 

In order to approximate fll.ip -f fLTl) we consider 

dt\ + div(<l>£(v)v (g) v) - div S(c, Dv) + Vp 

= -^^(div(Vc (g) Vc)) inVlx (0, T), (3.1) 

divv = 0, infix(0,T), (3.2) 

^tc + (^ev) ■ Vc = mA/i, infix(0,T), (3.3) 

= 0(c) - Ac. in X (0, T) (3.4) 

together with fll.5p -f lTT7|) . where \&eW = Po-(^£ * ^)\n, i^ei^) = £~'^ip{x/e), 
e > 0, is a usual smoothing kernel such that il){—x) = ip{x) for all x G M", w 
is extended by outside of Q, and Pa is the Helmholtz projection. Moreover, 
$s(s) = $(e|s|2) for all s G e > with some $ G Co"^(M) with $(0) = 1. 

The approximate system is formulated weakly as follows: For any ip G 
C^iQrY with divcp = and supp((^) CC n x [0,T), 

-/ w ■ dtLp d{x,t) + I (S(c,Dv) - V® v$e(v)) : Dcp c/(x,t) 

Jqt Jqt 

= - I ^£ (div(Vc(8) Vc)) ■ rf(x,t) + / WQ-ifi{Q)dx, 

holds and for any ip G C°°(n x [0,T)) with supp((^) CC H x [0,T), 
cdtip d{x,t) - / coV^(0)(ix+ / (^eV ■ Vc)^/' c?(x, t) 

= -/ mVjj^-Vij d{x,t) (3.5) 
n-Vc|af7 = (3.6) 

holds and 

= (y9(c) — Ac in Qt- 
For the following let Vp{n) = W^oinY n LP{n). 

THEOREM 3.1 Let ^ < q < oo, d = 2,3. For every < T < oo, 
Vq G L'^{Q), Co G i/^(i7) such that Cq{x) G [a, &] almost everywhere there is a 
weak solution (v,c, /i) of fl3.ip -f l3^ . fll.5p -f lL7|) such that 

V G l^J,(0, T; V;(fi)') n L^(0, T; V,m, 

c G C([0, T]; H\n)) n i7i(0, T; f^(-;(^])) H ^^(0, T; ^^^(j^))^ 

/i G L\0,T;H\n)) 
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where r = 6 if d = 3 and 1 < r < oo is arbitrary if d = 2. Moreover, for 
every <t <T 

lMt)\\l^n) + ErmMt)) + [ [ S(c,Dv) : Bwdxdr 
^ Jo Jn 

+ / m\V jj,]"^ dxdr = -\\vo\\l2m) + Emix{co) =■■ Eo (3.7) 



Jn 
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and 



\\c\\L2{0,T-W?{n)) + 110(c) II L2(o,T;L'-(f7)) < C{T,Eq) (3.8) 

for some C{T, Eq) > depending continuously on T, Eq. 
Proof: Let 

Xi := L'^{0,T-V,m n W;{0,T-V,{ny), Xo := L\0,T;Llm, 

and let Xi := [Xq, Xi)^iy Then Xi t- Xq by the Lemma of Aubin-Lions, 
cf. e.g. J.-L. Lions [15J or Simon tlTj, and since Xi ^ L°°(0, T; (fi)). 
Therefore 

Xi Xi Xo 

2 

due to [5l Theorem 3.8.1]. 

We define a mapping F : Xi — )■ Xi as follows: Given u G Xi , let c be the 

2 2 2 

solution of (EID-dM]) due to Theorem [O with v(x,t) = (^eu)(a;, t)x[o,T](^) 
and Co as in the assumptions. Then u c is continuous from the strong 
topology of Xq to the weak topology of 

Y = L\0, T- n H\Q, T- H^,]m 

as stated in Theorem 12. 1[ Therefore u i— )■ c is weakly continuous from Xi 
to Y. Moreover, Xo9ui— t-cGFR L°°(0, T; H^{Q)) is a bounded mapping 
and Y L'^{0,T; C^{Q)) by the Lemma of Aubin-Lions. Interpolation 

implies that Y 3 c ^ Vc e L^{0,T] L'^{n)) and y 9 c f-> Vc (g) Vc G 
L'^ {0,T] L^{Q)) are completely continuous mappings. Hence Xi 9 u h-> 

Vc ® Vc G lv'^'(0, T; L^(f2)) is completely continuous since Xi 9 u i— )• c G F 
is weakly continuous. 
Now let 

(f,^) = (f(u),(^) 
= / (Vc® Vc) : D(^s(^)ci(x,t) - / ($^(u)u ® u) : Dc^ c/(x, t) 



for all LP e L«(0,T; 1/^(1])) and let v = F(u) G Xi be the solution of the 
abstract evolution equation 

jvit) + A(t)v(t) = f (t) in L«'(0, T; (3.9) 

v(0) = vo mLlin), (3.10) 

where 

(A(t)(v(t)),(^)y,(n)',y,(Q)= / S(c(t),Dv(t)):D^da; 

for all G V^(f2). Then ^c(^) : K?(^) — ^ ^g(^)' is a strictly monotone, 
bounded, hemi-continuous, and coercive operator and existence of a unique 
solution follows from standard results on evolution equations for monotone 
operators, cf. e.g. [H^ Proposition 4.1]. Moreover, from standard estimates 
it follows that the mapping 

L«'(0, T; VgiQ)') X r 9 (f , c) ^ V G L^(0, T; Vg{Q)) H W^iO, T; VgiQ)') G X, 

is bounded. From this and the uniqueness of the solution one can derive that 
the latter mapping is also weakly continuous as follows: If 

(ffc,Cfc) -,^oo (f,c) in L^'{0,T;V,{Qy) x Y, 

then Cfc -^k^oo c in L^(0, T; C^(f2)) by compact embedding. Now, if is the 
solution of the evolution equation above with instead of f . Then (yk)km is 
bounded in Xi. Moreover, for any weakly convergent subsequence (vfcJjgN 
we have that 

A,/vfcJ - -^oo A(v) in L'>'{0,T- V,{Qy). 

Hence v^. ^j^oo v G Xi, where v is the unique solution of f l3.9p - fl3.10p . Be- 
cause of uniqueness of the solution, this holds true for any weakly convergent 
subsequence. Therefore ^fc^oo v G Xi. 

Moreover, it is easy to prove that Xi 9 u i— )■ f(u) G Xq is completely 
continuous. Altogether we obtain that the mapping Xi 9 u i— )■ v = -F(u) G 
Xi is weakly continuous. Therefore Xi 9 u h-> v = -F(u) G Xi is completely 
continuous. 

In order to apply the Leray-Schauder principle to F, cf. e.g. [TBI Chap- 
ter II, Lemma 3.1.1], it only remains to show that there is some R > such 
that 

AF(u) = u for some u G Xi , A G [0, 1] =^ ||u||xi < -R- 

2 2 
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Assume that AF(u) = u for some u G Xi, A G (0, 1]. (The case A = is 
trivial). Hence v = A"^u solves fl3.9p -( 1XTU]) with right-hand side /(u) as 
above. Thus taking the product of fl3.9p and v we conclude that 



^||v(T)||2+ / S(c,Dv):Dvf/(x,t) 

Jqt 



'Qt 

= ^llvolla - A"^(^'.u- Vu,u)q^ - (divVc® Vc,^',v)q^ 
= ^l|voil2 + A"'(/iVc,^',u)Q^ 
where we have used (11. 9p with e = 1. Combining this with (12. 5 p we obtain 

^||v(T)i|^+ / S(c,Dv):Dvt/(x,r) 
^ JQt 

+ -E,mix{c{T)) + -r i m|V/ip d{x, r) = -||vo||2 + t-^™x'(co) 
^ A Jq I A 



and therefore 



l|u|lL°°(0,T;L2(n)) + ll^llL'J(0,T;y9) 

= A^||v||^oo(o,T;L2) + -^'^ II U|| ^5(0^2.. y^) < Ci?(vo, Co). 

Because of (12. 6p . (12. 7p . there is some -R > such that 

||U|U, < M(||f(u)||^,,(o_^^^,)) < M'(||u|U^(o,T;L2)) < i?, 

where M, M' : [0, oo) — )■ [0, oo) are continuous and non- decreasing functions. 
Hence we can apply the Leray-Schauder principle to conclude the existence 
of a fixed point v = F{y\ v G Xi. Since the solution of (I3.9p -( 13TT0|) is in 
Xi, we even have v G Xi. Finally, the energy identity is proved by the same 
calculations as above with A = 1 and T replaced by t G (0,T). Finally it is 
easy to observe that (v, c, /i), is a weak solution in the sense above, where 
(c,ju) are determined by fl2TT|) - fl2l4l) with v(a;,t) = (\E'£u)(a;, t)x[o,T](i)- ■ 



Lemma 3.2 Let (£j)jeN a null sequence and (vj, Cj, ^j) be the solutions of 
(l3.ip - (P^ . (ll.5p - (ll.7p above with e replaced by Sj. Assume that vj — J-j-^oo v 
in L^(0,T; L'^{Q)). Then for a suitable subsequence 

Cj -,-^oo c m L\0,T; W^{n)), (3.11) 
c, ^j_,oo c m L^(0, T; Wl{n)), (3.12) 
/ij-j^oo/^ in L\0,T;H\Q)) (3.13) 
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where {c, fi) solve (11.31) . c\t=o = Cq, and n- V/i|aQ = in the sense that fll.l4p 
holds for anyip G C^iQx [0,T]) wi/isuppV^ CC Hx [0,T). Moreover, ffTij) 
/io/(is pointwise almost everywhere in Qt and n ■ Vc|aQ in dVt x (0,T) m ^/ie 
trace sense. 

Proof: Because of f l3.7p . fl3.8p . and fl3.4p . (cj)jer^ and are bounded in 

L^{0,T;W^{n)), L\0,T,H\n)), respectively. Hence (KTT\i and f l333|) hold 
for a suitable subsequence and some c G i^^(0, T; W^^(n)), yU G lv^(0, T; if ^(fi)). 
Moreover, (13. 3 p implies that (9tCj G L'^{0,T; H~^{fl)) is bounded. Hence 

c,-^,-^ooC mL\0,T-C\n)) 

by the Lemma of Aubin-Lions and W^{fl) C"'^(r2). Using that (cj)jggN C 

-Bf/C([0, T]; is bounded, a simple interpolation arguments yields (I3.12p . 

In order to prove (II. 4p . we use a monotonicity argument. To this end 
let 0o(s) := + as for all s G M and let A: ViA) C L\n) L\n) be 
defined A{c) = -Ac + 0o(c) for all c G V{A) with 

P(^) = {m G : 0oH e L2(fi),0[,(M)|VM|2 G L\n),dnu\9n = 0} 

Then 0o • (ct) 6) — M is monotone, ^ is a maximal monotone operator, and 
there is some C > such that 

MnHu) + \\Mc)\\LHn) < C (M(c) |U2(f^) + ||c|U2(f^) + l) (3.14) 

for all c G 'D{A) because of [l, Theorem 3.12.8], which is a variant of [U 
Theorem 4.3]. 

Therefore At: V{At) C ^^(fi x (0,T)) ^ L2(1] x (0,T)) defined by 
{ATc){t) = A{c(t)) for almost every t G (0,T) and c in 

^'(^t) = {uEL'^inx (0, T)) : u{t) G I?(^) for a.e. t G (0, T), 
^(m(-)) G L2(fi X (o,r)} 

is a monotone operator. Since ^ is maximal monotone, I + A: T>{A) — )■ 
L^(fi) is bijective. Combining this with (I3.14p . one easily obtains that also 
I + At: V{At) L'^{VL x (0,T)) is bijective, cf. e.g. [IS Lemma 1.3, 
Chapter IV]. Moreover, because of (13. 4p . 

At^Cj) = 4>o{cj) - Acj ^j^oo fJ'- ac in L^(fi x (0, T)). 

Furthermore, 

/ / ifJ-j — ctCj) Cj dx dt ^j-^oo / / {^i — ac)cdxdt 

Jo Jn Jo Jn 
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since Cj -^j^oo c and /ij ^j^oo in L'^{^ x (0,T)). Hence 

/ / A{cj)cj dx dt ^j^oo / / {^i — oic)cdx dt 
Jo Jn Jo Jn 

and [ini Proposition IV. 1.6] implies that c G P(.4t) and 

v4t(c) = /j- — ac, 

which is equivalent to (11.41) . 

Finally, fll.l4p follows easily by passing to the limit in (13. 5p and using the 
fact that 

for all V G L'^{0, T, L^(fi)). Furthermore n ■ Vc\dn = ^i^j-^oo n • Vcj |an = 
by the continuity of the trace operator. ■ 



4 Weak Solution 

To construct a weak solution of the above system, we use solutions of the 
approximate system (I3.ip - (I3.4I) together with (ll.5l) - (ll.6p . The existence of 
weak solutions of the approximate system follows from Theorem 13. 1[ In the 
following the solutions of the latter system are denoted by (ve,C£,/ie) for 
e > 0. Using the a priori estimates given by (13. 7p and (13. Sp . we can conclude 
for a suitable subsequence Ei — )-j^oo that 

Dv^. — )■ Dv weakly in L'^{Qt), 

n-\-2 

Vg. — )■ V weakly in L'^~^ (Qt), 
S(q„Dv,J^S weakly in /."'(Qt), 
® ^e^'^e, (v^J ^ H weakly in ^—{Qt). (4.1) 

Moreover, because of (13.81) . (13. 3p . and the Lemma of Aubin-Lions, it is easy 
to prove that 

Ce.^i^ooC mL\Q,T-C\Tl)) 

since W|(n) ^ C^{VL) compactly. Interpolation with the boundedness of 
Ce G L°^{0,T; H\n)) yields 

Vc,, ^i^oo Vc in L\Qt). (4.2) 
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Let K, e L'^{QtY'"^ be such that 

/ Ke : Dv? t)= I Vce ® Vce : 'D'^eiV') d{x, t) 
Jqt JQt 

-I Vc^Vc-.Bcp d{x,t) (4.3) 
for all cp G ^^(0, T; H^inY) and that K, G L'^{QtY and 

||K,||i2(Q^) < C||div(Vc® Vc) -^eCiiv(Vc,® VCe)||i2(o,T;/f-i) 

for some C > 0. We can assume that is pointwise a symmetric matrix. 
Then 

K,^-^0 strongly in L\QtY''\ 

due to (|i2D. 

Since q > J^, there exists some ctq > 1 such that > ctq > 1. Hence, 
due to fl4.ip we have for some -^i^oo 0, 

V strongly in L'^'"^{Qt) (4.4) 
and (g) V£^$£^(|ve|) V (g) V strongly in L"°{Qt). (4.5) 

We also have for i — )■ oo, 

V£^->v strongly in L''(0,T;L^(fi)), for all 1 < r < oo (4.6) 

by interpolation of (14. 4p with the boundedness of (ve)£g(o,i) G L°°(0, T; L^(r2)). 

Taking the limit of the weak form of the approximate system along the 
subsequence e^, we obtain the following limit equation: 



v-dt(p d{x,t)+ {S-v(S)v) -.Bcp d{x,t) (4.7) 
Jqt 

Vc (g) Vc : I)(p d{x, t)+ / vo • (fiiO) dx. 



for all ip G C^iQxY with div (p = 0. 

By subtracting the above equation from the weak form of the approximate 
equations, we have the following. 



/ {w,-^)-dt^ d{x,t)+ f (s(c„Dv,) -S) : Dcp rf(x,t) 
Jqt Jqt ^ ^ 

-- / (v, ®v,$e(v,)-v®v) :D(^c/(x,t)+ / K, : D(^rf(x, t). 

Jqt Jqt 
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Define := — v, then we can write this as 



J u^^ ■dt(pd{x,t) = y" : V(pd{x,t) 



(4i 



Qt Qt 

for any cp G C°°{QtY ^i^h div<^ = and supp(v3) CC VL x [0,T), where 
H, := Hi,i + H2,. with 

Hi,i := S(ce,,DveJ - S, 

H2,i := V,, ® v,,$,,(v,J - V ® V - K,, 

Then we have the following convergences for suitable Si — T-i^oo 



u^, weakly in L«(0,t; Vp{VL)), 

^ strongly in L^'"'{Qt), 

u,^ *-weakly in L°°(0, T; ^^(fi)), 

Hi,i^O weakly in L«'(Qt), 

Hi,2^0 strongly in L'"o(gT) 



(4.9) 
(4.10) 
(4.11) 
(4.12) 
(4.13) 



for some 1 < o"o < niin(g, g'). Let Jq be a time interval, Bq C R'^ be a ball 
such that Qo := h x Bq CC Qt- Let C e Co"^(|Qo) with xiq^ < C < Xiqq- 
Then we can apply Theorem 12.21 with K = S — S(c, Dv) to obtain 



lim sup 



((Hi,i + S - S(c,Dv)) : V(v,^ - v))CXoC dix,t) 



In other words 



lim sup 



1 ((S(c,„Dv,J-S(c,Dv)) :D(v,,-v))CXc.c^rf(x,t) 



Let G (0, 1). Then by Holder's inequality and < C < 1 we get 



lim sup 



1 ((S(c,,,Dv,J-S(c,Dv)) :D(v,,-v))'cXo,.c/(x,t) 
< c lim sup |Cj,fc|"'^ ^ < c2^*-"'^^^^^. 
This, the previous estimate and Holder's inequality give 



lim sup 



1 ((S(c,,,Dv,J-S(c,Dv)) :D(v,, -v))'crf(a:,t) 



<c2-(i-^ 
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For — 7- oo the right hand side converges to zero. Now, the monotonicity of 
S and ( > X|qo imphes that S(ce.,DveJ — )■ S(c, Dv) a.e. on ^Qq. Since Qo 
was an arbitrary space time cyhnder in Qt, we get S(ce;,Dv£.) — )■ S(c, Dv) 
a.e. on Qt- Since S(c£,,Dve.) — )■ S(c, Dv) weakly in L'^ (Qt), we get as 
desired S = S(c, Dv). This and (14. 7p prove that v solves 

-/ Y-dt(pd{x,t)+ (S(c,Dv) - V® v) : Dtp c/(x,t) (4.14) 
Jqt JQt 

= / K -.Bip d{x,t) + / vo-ifi{0)dx 

Jqt jo 

for any cp G C°°{QTy with divc^ = and supp((^) CC n x [0,T). 
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